Behaviour of relativistic black hole accretion sufficiently close to the horizon 
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Abstract 

This work introduces a novel formalism to investigate the role of the spin of astrophysical black 
holes in determining the behaviour of matter falling onto such accretors. Equations describing 
the general relativistic hydrodynamic accretion flow in the Kerr metric are formulated, and 
stationary solutions for such flow equations are provided. The accreting matter may become 
multi-transonic, allowing a stationary shock to form for certain initial boundary conditions. 
Such a shock determines the disc geometry and can drive strong outflows. The properties of 
matter extremely close to the event horizon are studied as a function of the Kerr parameter, 
leading to the possibility of detecting a new spectral signature of black hole spin. 



1 Flow Structure 
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The energy momentum tensor Q^'^ has been formulated (the general most form of 9^*^ is available 
in Novikov & Thorne 1973) in the Boyer-Lindquist co-ordinates, and its covariant derivative has 
been evaluated to obtain the general relativistic Euler equation and the equation of continuity, 
in the form of a set of spatio-temporal first order differential equations. Axi-symmetry has been 
adopted for the flow geometry as described in figure 1, where, following Abramowicz, Lanza & 
Percival (1997), the disk height h{r) is calculated as: 
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Two first integrals of motion for the system defined along streamlines, the conserved specific flow 
energy, £ and the baryonic load rate, M, along with the corresponding entropy accretion rate, H, 
are calculated as: 
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Differential solutions of the first integrals of motion provides the three velocity gradient (of the 
flow) as a first order autonomous dynamical system: 
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In Eqns. (1-4), r, u, Cs,p, A, 7, a are the radial distance on the equatorial plane, three (flow) veloc- 
ity, sound speed, density, flow specific angular momentum, flow adiabatic index and the black 
hole spin (the Kerr parameter) respectively. A and A are functions of r and a, Vf is function of 
u, r, A, a, A, A, and other terms (especially in (4)) are functions of above quantities and of various 
metric elements as well as their space derivatives. Further details about the above mentioned 
equations are available in Barai, Das & Wiita (2004). 



2 Transonicity, Shock and Outflow Generation 

The fixed point solution obtained from (4) provides the number of critical points the flow can pass 
through, and thus provides the number of transonic flips. 

Region of multi-transonicity for both accretion (marked A) and wind (marked W) in the pa- 
rameter space are shown in figure 2. The labels I and O indicate regions with lone inner and outer 
critical points, respectively. The plot has been generated for a = 0.3 and 7 = 1.33. 

Multi-transonic (actually, bi-modal, since the middle sonic point does not support a steady 
solution that passes through it) accretion flows may encounter a stationary shock. The general 
relativistic shock conditions are formulated: 

[[puT^]] = 0, [[Qti^v"]] = [[{P + e)vtuT^]] = 0, [[SJ^.r/'^r?^] = [[{? + e^Tl +p]]=0, (5) 

where T^ is the Lorentz factor, and rj^ is the normal to the hypersurface S of discontinuity. The 
expression for the shock invariant Sh is computed: 
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The shock strength and the entropy enhancement ratio (at the shock) are defined as the ratios of 
Mach numbers (M_/M_|_) and the entropy accretion rate (S_|_/iE_) respectively, where — and -h 
indicate the pre- and post-shock values of any flow variable. 

Simultaneous solution of Eqns. (1-6) provides the phase portrait of the shocked multi-transonic 
accretion, as shown in figure 3. Transonic flow (marked by 'A') through the outer sonic point en- 
counters a stable stationary Rankine-Hugoniot shock transition marked by SI (S2 being the other 
formal shock location, which is unstable and hence is not considered), produces post shock sub- 
sonic flow, which again becomes supersonic (segment marked by 'a') after passing through the 
inner sonic point; and finally plunges through the event horizon. As a consequence of the shock 
formation, the post shock flow has higher temperature, density, pressure and residence time com- 
pared to its pre-shock coimterpart. This favours the generation of an optically thick halo (yellow 
coloured elliptic patches in figure 4, where the disc structure has been obtained by solving (1-6)) in 
the post shock region as a result of disc evaporation, and initiates the production of thermally and 
centrifugally driven cosmic outflows. 

Figure 5a shows the variation of shock location rgh as a function of black hole spin. Figure 5b 
shows that the stronger shocks (higher value of (M_/M_|_)) are produced closer to the black hole, 
and is a result of maximum entropy production (measured by (H+/H_)) due to the shock. Pre- 
and post-shock ratios of temperature, pressure, density and velocity (some of the normalized by 
appropriate numerical factors so as to fit in the same figure) as a function of black hole spin are 
plotted in figure 5c. 



3 Quasi-terminal Values 

Flow variables (calculated along the solution 'a') at a very close proximity = r+ -I- 5 (r+ = 
1 -I- -\/(l — a^), and 5 = O.OOlrp) of the event horizon are termed as quasi-terminal values, and are 
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distinguished with a subscript S. The variation of the (normalized) quasi-terminal temperature 
(T^), density (ps), pressure (p^), as a function of the black hole spin, are shown in the figure 6. The 
calculations are performed for a 3x IO^Mq black hole accreting at a rate of 4.29x lO~^M0Yr~^. 

Shock formation is naturally also possible for other values of the black hole mass and the accretion 
rate, as well as for values of the black hole spin higher than those shown here in Fig. 6. Accretion 
onto the rotating black hole with the Kerr parameter as high as 0.95 has been observed to form 
shocks; the details of such work will be reported elsewhere (Barai, Chakraborty, Das & Wiita, in 
preparation). 

4 Conclusion 

• For the first time, the shocked relativistic accretion in the Kerr metric has been studied at such 
a close proximity to the black hole event horizon. 

• The effects of the black hole spin on the dynamical and thermodynamical properties of the 
accreting material have been made explicit. 

• Accurate estimates of the flow temperature, velocity and density profile has been made. This 
leads to a possible formalism for the study of the spectral signature of the black hole spin 
(work in progress). 
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Figure 1: The yellow circular patch with BH written inside represents the black hole and the pink 
dashed boundary mimics the event horizon. The wedge shaped dark green lines represents the 
envelop of the accretion disc. The light green line centrally flanked by the two dark green disk 
boundaries, is the equatorial plane, on which all of the dynamical quantities (e.g., the advective 
three velocity u) are assumed to be confined. Any thermodynamic quantity (e.g., the flow density) 
is averaged over the local disc height h(r). 
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Figure 2 
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Figure 3 



SHOCKED INNER PART OF ACCRETION DISC AROUND A KB<R BLACK «OLE 



300 



200 



100 



PRE-SHOCK POST- - ' POST- 

SHOCK , . ' ' ■ - - SHOCK 



-100 



-200 



-300 



PRE.-SHOCK 



-150 



-100 



-SO 



so 



100 



ISO 



Figure 4: Pre- and post-shock disc geometry with thermally driven optically thick halo. 
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Figure 5a 
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Figure 5b 
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Figure 5c 
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Figure 6a 
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Figure 6b 
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Figure 6c 
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